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The so-called black hole shadow is a dark region which is expected to appear in a 
fine image of optical observation of black holes. It is essentially an absorption cross 
section of black hole, and the boundary of shadow is determined by unstable circular 
orbits of photons (UCOP). If there exists a compact object possessing UCOP but no 
black hole horizon, it can provide us with the same shadow image with black holes, 
and a detection of shadow image cannot be a direct evidence of black hole existence. 
Then, this paper examine whether or not such compact objects can exist under some 
suitable conditions. We investigate thoroughly the static spherical polytropic ball of 
perfect fluid with single polytrope index, and then investigate a representative example 
of the piecewise polytropic ball. Our result is that the spherical polytropic ball which 
we have investigated cannot possess UCOP, if the sound speed at center is subluminal 
(slower-than-light). This means that, if the polytrope treated in this paper is a good 
model of stellar matter in compact objects, the detection of shadow image is regarded 
as a good evidence of black hole existence. As a by-product, we have found the upper 
bound of the mass-to-radius radio of polytropic ball with single index, < 0.281, 

under the subluminal-sound-speed condition. 

Subject Index E31, E30, E01 


1. Introduction and our aim 

In recent years, the resolution of image by the very-long-baseline-interferometer (VLBI) radio 
observation is approaching to the visible angular size of SgrA*, about 10 micro-arcsecond (a 
black hole candidate of 4 x 10 6 M 0 at the center of our galaxy, 8 kpc from the earth), which 
is the largest visible angular size in known black hole candidates [1, 2]. The so-called black 
hole shadow is expected to be resolved by such fine observation near future (see [3, 4] and 
references therein). It seems to be currently a common understanding that seeing the black 
hole shadow is believing the existence of black hole horizon. 

However, this common understanding has not been confirmed in general relativity as fol¬ 
lows: Remember that the black hole shadow is a dark region appearing in an optical image 
of black holes, on which some photons would be detected if the black hole did not exist. 
Therefore, the shadow is essentially an absorption cross section of black hole. However, it 
should be emphasized that photons on the edge of shadow have been circulating around 
black hole before coming to the observer. The innermost circular orbit of those photons is 
not a great circle on black hole horizon, but an unstable circular orbit of photons (UCOP). 
That is to say, the boundary of shadow is determined not by the black hole horizon, but 
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by the UCOP. This indicates that the direct origin of shadow is the UCOP, not the black 
hole horizon. Hence, although we can conclude the existence of UCOP once a shadow is 
observed, however, we cannot conclude immediately the existence of black hole horizon even 
if a shadow is clearly observed. 

Here, let us assign a term, black hole mimicker, to a compact object possessing UCOP 
but no black hole horizon. If there exists a black hole mimicker, it can provide us with 
the same shadow image with black holes in optical observation, and a detection of shadow 
image cannot be a direct evidence of black hole existence. Therefore, we are interested in an 
existence/non-existence condition of black hole mimickers. 

Some exotic candidates of black hole mimicker have been proposed, such as gravastars and 
boson stars. Under the assumption that any black hole mimicker emits a thermal radiation 
from its surface, the current observational data of SgrA* (the mass accretion rate and the 
observed flux) exclude the possibility that SgrA* may be those exotic black hole mimick¬ 
ers [5]. Further, the current observational data of some black hole candidates indicate that 
the gravastar cannot exist in nature, if any gravastar emits always a thermal radiation from 
its surface [6]. For the case that the gravastar does not emit radiations from its surface, the 
general feature of gravastar shadow, which enables us to distinguish gravastars from black 
holes in the shadow image, has already been examined [7]. The investigation on those exotic 
models may be interesting. However, we focus on a rather usual model in this paper. 

Consider a static spherical ball of perfect fluid matter in the framework of general rela¬ 
tivity, which connects to Schwarzschild geometry at its surface. A fluid ball, which does not 
possess the black hole horizon, becomes a black hole mimicker if it possesses one of following 
properties: 

(A) The ball is so compact that there appears a UCOP in the outside Schwarzschild 
geometry. 

(B) The ball is not so compact as case (A), but a UCOP appears inside the ball. 

In case (A), if the surface of fluid ball neither emit nor reflect any radiation, this ball can 
provide us with the same shadow image with a black hole. In case (B), if the fluid outside 
UCOP is completely transparent and if the fluid inside UCOP is not transparent, this ball 
can provide us with the same shadow image with a black hole. 

Concerning the case (A), the mass-to-radius radio 3M*/I7* of fluid ball is the key quantity, 
where M* and R* are respectively the total mass and surface radius of the ball measured 
in the dimension of length. If this ratio is less than unity (3M*/.R* < 1), then no UCOP 
appears outside the fluid ball because the radius of UCOP in Schwarzschild geometry is 
3M*. Note that, in order to avoid the gravitational collapse, an inequality, 3M*/i?* < 3/2, 
must hold. Further, by adding some reasonable conditions to the fluid ball, the upper bound 
of the ratio, 3M*/i?* < U , should decrease, 


3M» 

R* 


<U <-. 
2 


( 1 ) 


Hence, the problem in case (A) is whether or not the upper bound U becomes less than 
unity under some reasonable conditions of the fluid ball. 

The upper bound U has already been estimated for some situations. For the fluid ball with 
any equation of states satisfying three conditions, non-increasing energy density in outward 
direction, barotropic form of equation of states, and subluminal (slower-than-light) sound 
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speed, an upper bound has been obtained by our previous work [8], U ~ 1.0909209. And, 
for a core-envelope model of neutron star [9, 10], the other value of upper bound has been 
obtained, U — 1.018, where the equation of states in the envelope region is determined by 
nuclear matter physics (see section 2 of [10] for details), while the equation of states in the 
core region is treaded by the same method used in our previous work [8]. These estimations 
give still a bound greater than unity, U > 1. Further investigation is necessary to consider 
what kind of equation of states results in U < 1. 

On the other hand, concerning the case (B), there is no existing work analyzing UCOP 
inside a fluid ball as far as we know. We need to formulate the criterion for judging the 
existence/non-existence of UCOP, and then apply the criterion to our model of fluid ball. 

In this paper, we regard the polytrope as a representative example of the barotropic 
equation of states. The model investigated here is the static spherical ball of perfect fluid with 
the polytropic equation of states. Our investigation focuses mainly on the simple polytrope 
model whose polytrope index is fixed at one value for all region inside the fluid ball. After 
the thorough investigation of the simple model, an example of the piecewise polytropic 
fluid ball [11] is investigated in the frame work of core-envelope model, where the value of 
polytrope index in the core region differs from that in the envelope region. Although the 
thorough study of piecewise polytrope model is left for next tasks, we can find a good insight 
into the core-envelope model of polytropic fluid balls. 

For both polytrope models, we examine whether or not the cases (A) and (B) are possible 
for the polytropic fluid ball. Our result is that there cannot exist UCOP in neither outside 
nor inside of the polytropic fluid balls, if the sound speed at center is subluminal. This 
implies that, if the polytrope investigated in this paper is a good model of stellar matter in 
compact objects, a detection of shadow image is regarded as a good evidence of existence of 
black hole. 1 

In section 2, a set-up of our analysis is described. Sections 3 and 4 are devoted to analyses 
of, respectively, cases (A) and (B) for the simple polytrope model. In section 5, the analyses 
in previous sections are extended to an example of the core-envelope piecewise polytrope 
model. Section 6 is for summary and discussions. 

2. Static spherical polytropic ball 

We consider static spherical ball made of polytropic perfect fluid. The metric of this 
spacetime is given by a line element, 

dr- 2 

ds 2 = g^ u dx^ dx v = -e 2$(r) c 2 df 2 + -—— 2 + r 2 (dd 2 + sin 2 6d<p 2 ) , (2) 

1 — 2Gm(r)/(c^r) 

where (t, r, 9, tp) is spherical poler coordinates, 4>(r) gives the lapse function, and m{r ) is the 
mass of fluid contained in the spherical region of radius r. The energy-momentum tensor of 
perfect fluid is T^ v = [ a(r)c 2 + p{r)]u^u v + p(r)g^, where u = e~®d ct is the four-velocity 
of static fluid, and cr(r) and p{r ) are respectively the mass density and pressure of fluid. 


1 There have been some attempts to study a nonlinear instability of black hole mimickers (see [12] 
and references therein). This is an interesting approach. However, in those works [12], while the 
argument for nonlinear instability of black hole mimicker is conjectured by combination of linear 
analyses, but no definite proof of nonlinear instability has been obtained. 
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A condition, m(0) = 0, should hold due to the regularity of spacetime at center. This 
implies a finite mass density at center, a c = ct( 0 ) 7 ^ 00 , where the suffix c denotes the value 
at center. We normalize all quantities by a c , 


R VG^c 

P := - r 


m) ■■= 


cr(r) 


M(R) := 




(7c 


m(r ) , P(P) '■= 


p(r) 


(3) 


C (Jq C CT qC 

These are dimension-less. The lapse function, < J>(r) := <h(P), does not need normalization 
because $ is already dimension-less by definition ( 2 ). 

The barotropic equation of states is generally expressed as P = P(S). We adopt the 
polytrope as a representative form of the barotropic matter, 


P(E) = PE 1+1/n , 


(4a) 


where K and n are positive constants, and n is called the polytrope index. By normaliza¬ 
tion (3), the mass density at center is unity, E c = 1. Therefore, the coefficient K is equal to 
the pressure at center in our normalization (3), 

Pc = K (= P(E = 1)). (4b) 

The form (4) is the simple polytrope whose index, n, is fixed at one value for all region 
inside the fluid ball. We are going to extend the simple form (4) to the piecewise polytrope 
in section 5, but in this section we focus on the simple form (4). 

The surface of fluid ball is defined by the zero pressure, P* = 0, where the suffix * denotes 
the value at surface. Therefore, the mass density at surface is zero due to polytropic equation 
of states (4), E* = 0. The normalized mass density takes values in the interval, 

(£* =)0<S< 1(=E C ). (5) 


Regarding E as an independent variable, this finite interval of E seems to be useful for 
numerical analysis. Hence, we regard all variables as functions of E, 

P = P( E) , M = M( E) , $ = $(E), ( 6 ) 

and P(E) is already expressed as a function of E in (4). The surface radius P* and total 
mass M* of polytropic ball are determined by 

P* = p(E = 0) , M* = M(E = 0). (7) 


The sound speed V in polytropic ball is given by 


2 = dP(S) 
dE 


= P c (l + i)E 1 / n . 


( 8 ) 


This sound speed is normalized by light speed. Obviously, this V decreases from center 
(E c = 1) to surface (E* = 0). The highest sound speed is given at the center, 

1 


= P 


K)- 


(9) 


In this paper, we assume a subluminal condition of sound speed, 


V c < 1. (10) 

Here it may be fair to note that a possibility of omitting the subluminal condition (10) 
has been discussed for the nuclear matters of neutron stars [9]. The omission of condi¬ 
tion ( 10 ) may be possible if the nuclear matter possesses the properties so that the matter 
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temperature is zero and the sound wave is dispersed and dumped (absorbed) by matters 
instantaneously. Even though such properties are good approximate ones of nuclear matters 
in neutron stars, those properties may not be exact and correct for all region inside neu¬ 
tron stars. Furthermore, our analysis is not restricted to neutron stars but is designed to 
include any matter described by polytropic equation of stats under some conditions which 
are reasonable from the viewpoint of general relativity. Hence, let us require the subluminal- 
sound-speed condition (10) which we regard as a rigorous general relativistic property of 
matters. 

The outside region of polytropic ball, P > P*, is described by the Schwarzschild geom¬ 
etry of mass M*. The inside region P < P* is determined by the Einstein equation and 
conservation law T^ LU - U = 0, which are reduced to the Tolman-Oppenheimer-Volkoff (TOV) 
equations, 


d M 
dS 
dP 
dE 
dd> 
dE 


(S + P) (M + 4vrP 3 P) d R 


P(P - 2 M) 
1 d P 


dE 


E + PdE 


(11a) 

(lib) 

(11c) 


Two functions P(E) and M(E) are obtained by solving (11a) and (lib) under the equation 
of states (4). Substituting those solutions into (11c), <I>(E) is obtained. Those solutions of 
TOV equations depend on two parameters, V c and n, for the simple polytrope (4). 

Under the set-up given above, our aim is to analyze the problem, whether or not the 
properties (A) and (B), which are described in section 1, hold for polytropic ball under 
the subluminal-sound-speed condition (10). In following analyses, TOV equations (11a) 
and (lib) are solved numerically. A technical remark for numerical calculation is summa¬ 
rized in the appendix A, which is applied to sections 3, 4 and 5. All of our numerical analyses 
are performed with Mathematica ver.10. 


3. Problem A: Can a UCOP appear outside the simple polytropic ball? 

The problem in this section is whether or not an inequality, P* < 3M*, holds for the simple 
polytropic ball of equation of states (4) under the condition (10). If P* < 3M*, then a UCOP 
appears outside a polytropic ball. Otherwise, if 3M* < P*, then a UCOP does not appear 
outside the ball. Our strategy is as follows: 

Al: Solve numerically TOV equations (11a) and (lib) for given values of parameters (14, n), 
and calculate the mass-to-radius ratio, 3M*/P* . 

A2: Iterate the step Al with varying parameters (14, n), so as to obtain the ratio, 3M*/P*, 
as a function of parameters (14, n). 

A3: Find the maximum value of 3M*/P* as a function of (14, n). If the maximum is less 
than unity, we conclude that the inequality, 3M* < P*, holds for all values of (14, n), 
and no UCOP appears outside the simple polytropic ball of equation of states (4) under 
the subluminal-sound-speed condition (10). 

In the Newton gravity, the total mass and surface radius of the simple polytropic ball 
are finite in the interval of index, 0 < n < 5, but diverge in the interval, 5 < n, for any 
value of 14 > 0. However, in the Einstein gravity, the thorough numerical analysis of simple 
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Fig. 1 Numerical solution of TOV equations with the simple polytropic equation of 
states (4) at (V c ,n) = (0.6,6). Although M* and R * are infinity, the mass-to-radius ratio 
M*/becomes zero. 


polytropic ball by Nilsson and Uggla [13] have revealed a complicated behavior of M* and 
R* in the half-infinite interval of central sound speed, 0 < V c : 

o In the interval of polytrope index, 0 < n < 3.339, both of M* and I?* are finite, 
o In the interval, 3.339 < n < 5, a complicated behavior is found. 

• Both of Mif and R* are finite for almost of all values of (V c ,n) in the present parameter 
region. 

• However, both of M* and ii* diverge at some discrete points (Vq^ , n®) in V c -n plane, 
where i = 1,2,--- ,N ^ iv . The number of such divergence points, A r c j iv (= finite or 
countable infinity), cannot be read from Nilsson-Uggla [13]. 

Note that, although the existence of some divergence points (Vc l \n^) in V c -n plane 
has been definitely confirmed, the accurate positions of them have not been specified. 2 
o In the interval, 5 < n, both of M* and R * are infinity. 

Note that the mass-to-radius ratio has not been analyzed in Nilsson-Uggla [13]. The analysis 
of the ratio is our task. 

From the above behavior of M* and R * found by Nilsson-Uggla, the physically interesting 
region of parameters are 

0 < U c < 1 , 0 < n < 5 , (12) 


2 Two examples of such divergence points are (14, n) = (I4^\ 3.357), (u) 2 \ 4.414), where the values 
U' 1 ^ and Vc 2 ' 1 cannot be read from Nilsson-Uggla [13]. 
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contours of 


3M* 



Vc 


Fig. 2 Contours of 3on V c -n plane for the simple polytropic ball. Numerical calcula¬ 
tion is performed in the region, 0.05 < V c < 1.5 and 0.01 < n < 5.5. Obviously, the inequality, 
3< 1, holds under the subluminal-sound-speed condition, V c < 1. For small values of 
14 and n, numerical errors become manifest. 


where the interval of 14 denotes the subluminal-sound-speed condition (10). It is enough for 
our aim to calculate the ratio, 3in this parameter region. However, the parameter 
points (14where M* and i?* diverge in the region (12), may not be included in 
the grid points of numerical analysis (see the step A2 of our strategy). In order to guess the 
behavior of 3at those points (14^,n^), we observe the solutions of TOV equations in 
the interval, 5 < n, where M* and R* diverge as well. Figure 1 is an example with 14 = 0.6 
and n = 6. This figure shows that, although the mass M(E) and radius i?(E) diverge as 
the surface (E = 0) is approached, the mass-to-radius ratio M(E)/i?(E) converges to zero. 
The same behavior is observed for the other values of (14, n) in the interval, 5 < n. Hence, 
it is expected, even at the parameter points (14^, n^) where M* and R * diverge in the 
region (12), the mass-to-radius ratio converges to zero. 

With the help of above discussion, we can safely carry out our strategy of numerical analy¬ 
sis, composed of steps Al, A2 and A3. The result is shown in figure 2, in which the contours 
of 3are plotted. Although our main interest is in the parameter region (12), we have 
calculated for a bit larger region, 0.05 < 14 < 1.5 and 0.01 < n < 5.5. 3 It is obvious that the 


3 In figure 2, the ratio 3is small enough in the interval, 5 < n < 5.5. This is consistent with 
the figure 1 and discussion after (12). Further, we have found numerical implications, although details 
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If UCOP exits, then ... 



Fig. 3 A schematic graph of U e s(R) if a UCOP exists inside the polytropic ball. 


maximum value of 3in the region (12) appears on the vertical line at V c = 1 in figure 2. 
This maximum takes the value between 3M*/R* = 0.8 and 0.85. (A more precise value is 
calculated in section 6.) Hence, we can conclude that the inequality, 3 M*/R* < 1, holds 
in the physically interesting region (12). No UCOP appears in the outside Schwarzschild 
geometry under the subluminal-sound-speed condition (10) for the simple polytropic ball of 
equation of states (4). 

4. Problem B: Can a UCOP appear inside the simple polytropic ball? 

The problem in this section is whether a UCOP can exist inside the simple polytropic ball 
of equation of states (4) under the subluminal-sound-speed condition (10). For the first, the 
conditions for the existence of UCOP in static spherical spacetimes are summarized. Next, 
those conditions are applied to the simple polytropic ball, and we show that no UCOP can 
exit inside the simple polytropic ball. 


4-1- UCOP in static spherical spacetimes 

The feature of UCOP in the static spherical spacetime of metric (2) is determined by the 
null geodesic equation. Denoting the affine parameter and radial coordinate of null geodesic 
by, respectively, A and i? n uii(A) under the normalization (3), the radial component of null 
geodesic equation on spacetime (2) is 

(4^) 2 +^!/dT(^„ 1|) = 0, (13) 

where U e s is the effective potential given by 


U eS (R) 


\— 

-R 2 


exp(—24>(i?)) 



2 M(R)\ 

R~ ) 



(14) 


where M(R) and Q(R) are regarded as functions of R given by solving TOV equations, b is 
an impact parameter, and l and uj are respectively the orbital angular momentum and the 
frequency of photon measured at infinity. 


are not shown here, that the points (V c ^,nW), where M* and 1?* diverge in the region (12), form a 
line along the valley shown in figure 2, and this line seems to approach (V c , n) = (0, 5). 
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A photon propagating on a UCOP remains at a constant radius (dP nu n/dA = 0), but it is 
unstable. Hence, the radius of UCOP, R u , is determined by following conditions, 

U eff( R u) = 0 , ^-(R a ) = 0 , ^(i4)<0. (15) 


This implies that, if a UCOP exists inside the polytropic ball, the top of potential barrier 
touches below the zero level at R u as shown in figure 3. 


By substituting (14) into (15), we obtain 


= 0 

(16a) 

1 dP d 2 P 

(16b) 

P 2 

b 2 = P, 

F(Ru) ’ 

(16c) 

where 


F(R) :=exp(2cl»(P)) (= -g 00 ). 

(17) 


The radius of UCOP, R u , is determined by (16a) and (16b), and then the impact parameter 
of null geodesic circulating on UCOP forever is obtained by (16c). Therefore, the existence 
condition of UCOP consists of two parts; an algebraic equation (16a) and an inequality (16b), 
which do not include the impact parameter. 


4-2. Non-existence of UCOP inside the simple polytropic ball 

In order to apply the existence conditions of UCOP (16a) and (16b) to simple polytropic 
balls, we need a concrete functional form of <b(E) = (1/2) In P(E). Substituting the equation 
of states (4) into a TOY equation (11c), 


dd>(S) _ (i , 1 A PcS 1 /”- 1 _ ^ dlna + PcE 1 /") 

ds v + Ji + p c sv- ln+ j dS 


(18) 


The integration constant for this equation is determined by the junction condition of metric 
at the surface of polytropic ball, exp(24>*) = 1 — 2M*/P*. We obtain 


P(£) = e 2<I>(s) 


-p/i I ^Mr 2(n+1) 

(1 + P c S 1 /n)2(n+l) £ ) 


(19) 


where P* = 1 — 2M*/P* . Regarding S as a function of R, which is given by solving TOV 
equations (11a) and (lib), we obtain P(P) as a function of R. 

Prom (19), first and second differentials of F(R) are calculated, 


d F(R) M + 4ttR 3 P 

d R ^ R{R - 2 M) [ ' 


(> 0 ) 


d 2 F(R) _ 1 4 M-R 3 / S + P U + 3 P \ 

d R 2 R R-2M 71 \R-2M M + AttR 3 p) 


d F{R) 
d R 


( 20 ) 


where TOV equations (11a), (lib) and (11c) are used. Substituting these differentials into 
the existence conditions of UCOP (16a) and (16b), we obtain 


^ i | r=r u ~ 0 j C2 1 n=n n > 0, 


( 21 ) 


9/19 





















n = 0.7 


n = 1 


n = 3 






Fig. 4 Plots of curves Ci(£, V c ,n) = 0 and C 2 {Y>, , ri) = 0 in T,-V c plane for some fixed 

n. The regions, C 2 > 0 and C 2 < 0, are also denoted. It is recognized that the existence 
conditions of UCOP (21) are not satisfied for the single polytropic balls. 


where 


Ci = R - AnR 3 P - 3 M 


C 2 = 1- 


4 M — R g/S + P £ + 3 P \ 
R-2M \R-2M M + 4tvR 3 PJ 


( 22 ) 


If there does not exist R u which satisfies the conditions (21) for any value of parameters 
(Vc,n) in the physically interesting region (12), then it is concluded that no UCOP can 
appear inside the polytropic balls. 

Note that the value of quantities C\ and C 2 are calculated by substituting the solutions 
of TOV equations (11a) and (lib). The solutions of TOV equations are functions of £ and 
depend on parameters (V c ,n). Hence, in our analysis, C\ and C 2 can be obtained numeri¬ 
cally as functions of three arguments, Ci(£, V c ,n ) and C 2 (Y>, V c ,n). Then, in order to check 
whether or not there exists R u satisfying (21), our strategy is as follows: 

Bl: Solve numerically TOV equations (11a) and (lib) for given value of parameters (V c ,n), 
and iterate this numerical calculation with varying V c and fixing n at a given value. 
This iteration produces Ci(£,V c ,n) and C 2 (P,V c ,n) as functions of (£,U C ) for the 
given value of n. 


10/19 



















B2: Plot two curves, C\ = 0 and C *2 = 0, and identify two regions, C 2 > 0 and C 2 < 0, in 
Ti-V c plane for the given n. If the curve C\ = 0 does not intersect with the region C 2 > 0, 
then it is concluded that no UCOP exists inside the simple polytropic ball at the given 
value of n. 

B3: Iterate the steps B1 and B2 with varying n, and check whether or not the intersection 
of C± = 0 with C 2 > 0 exists at each value of n. If the intersection does not appear 
for any value of n, then we conclude that a UCOP can never appear inside the simple 
polytropic ball of equation of states (4). 

The numerical result of this strategy is shown in figure 4. We find that the curve C\ = 0 
remains in the region C 2 < 0 for all values of n = 0.7, 1, 3,4, 5, 6 . The same feature is 
found for the other values of n as far as we have calculated. Hence, no UCOP appears inside 
the polytropic ball. Furthermore, the non-existence of UCOP inside polytropic balls seems 
to hold for not only the physically interesting parameter region (12) but also all parameter 
region, 0 < n and 0 < V c . 


5. Extension to the core-envelope piecewise polytrope models 

In this section, the analysis of simple polytrope model performed so far are extended to 
a representative model of the core-envelope piecewise polytropic fluid balls. For the first, 
the formulation of our model is described in section 5.1. Then, in sections 5.2 and 5.3, we 
investigate whether or not the properties (A) and (B), which are described in section 1, hold 
for our model of the piecewise polytropic balls. Our analysis will indicate the same result 
with the simple polytropic ball that no UCOP appears inside and outside the representative 
model of the core-envelope piecewise polytropic balls. 


5.1. Our model of the core-envelope piecewise polytropic perfect fluid ball 
Dividing the interval of normalized mass density (0 < £ < 1) into some sub-intervals, the 
piecewise polytropic equation of states is defined as the one whose polytrope index takes 
a different value at each sub-interval [11]. The core-envelope type of it is the polytrope 
composed of two sub-intervals (see figure 5), 


P(£) 


P c £ 1+1 /"<- in the core region, £b < £ < 1 
iF env £ 1 + 1 /"- in the envelope region, 0 < £ < £b, 


(23) 


where £b is the mass density at the boundary between core and envelope regions, the 
parameters n cor and n env are the polytrope indices of respectively the core and envelope 
regions, and P c is the pressure at center. We require that the pressure is continuous at the 
boundary between two regions, P(£b — 0 ) = -P(£b + 0 ), in order to maintain the mechani¬ 
cal balance at the boundary. By this junction condition, the coefficient K cnv is given by the 
other parameters, 

K env = P c £V"--i/"~ . (24) 

This implies that dP/d£ becomes discontinuous at £b as shown in figure 5. And, through 
the TOY equation (lib), di?/d£ becomes discontinuous at £b- 
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Fig. 5 Left panel is a schematic graph of the core-envelope piecewise polytropic equation 
of states (23). Right panel is a schematic image of the core-envelope model of fluid ball. 


If we would restrict our discussion to neutron stars, the nuclear matter near the surface 
might be in a solid state, and the matter around the core might be in a fluid state [9]. 4 
Such inner structure of neutron star may be approximately modeled by applying the above 
equation of states (23), and the phase transition density Eb is determined by some reasonable 
nuclear matter physics. However, our discussion is not only for neutron stars but also for 
general spherical static balls of equation of states (23) under the subluminal-sound-speed 
condition. Therefore, we regard Eb as a free parameter in our analysis. 

The sound speed V is given by V 2 = d-P/dE. Then, there are two local maximums of V at 
the center (E = 1) and the boundary approached from the envelope (E —>• Eb — 0). Hence, 
the subluminal-sound-speed condition is given by 


Vr :=Pc 1 + 


n c 


< 1 


K := I<r 


1 + 


n e 


.1 /n.„v^pb + 1 


fl + —) 


1/n cor / -i 

J b — 4 ' 


(25) 


where the relation (24) is used in the second equality for Vb. 

Here let us note that, according to our previous paper [8] which analysed the mass-to-radius 
ratio 3M*/P* of generic perfect fluid balls of any equation of states under the subluminal- 
sound-speed condition, the ratio 3M*/P* tends to decrease as the sound speed increases 
inside the ball. Therefore, in search of the upper bound of 3M*/P* for the core-envelope 
piecewise polytropic balls, it seems to be reasonable to set the sound speed at E = 1 and Eb 
being the light speed, 

V c = l , Vb = l. (26) 


The polytropic ball under this extreme-sound-speed condition is the model we focus on in 
the following analyses. 


4 When going down from the surface to the center of neutron star, there may occur a phase 
transition from a solid state to a fluid state at some radius f?b = R( Eb) due to the increase of matter 
density in the inward direction. 
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The condition (26) and the relation (24) give the relations, 


-fc(^cor) 
^env ( ^ b 5 ^cor ) 
-Kenv (^b> ^ cor) 


Tin 


^cor T 1 


'cor 

^cor 


+ A 

- 


i -i 


- 1 


n , 


Tlcor 

n cor ^l/n cor -l/n e 

+ I b 

:or i 


(27a) 

(27b) 

(27c) 


Through these relations, we find that the core-envelope piecewise polytrope (23) depends on 
a variable X and two free parameters (Xb,n cor ). Actually, our equation of states (23) is the 
function of three arguments, T^X, X b , n cor ), under the extreme-sound-speed condition (26). 
The parameter Xb takes, obviously, a value in the interval, 


0 < X b < 1. 


(28a) 


On the other hand, the physically interesting interval of the piecewise polytrope index, n cor , is 
not obvious. Remember that, for the simple polytrope (4), the physically interesting interval 
of polytrope index (0 < n < 5) has been determined by Nilsson-Uggla’s thorough numerical 
analysis of the simple polytrope [13]. Therefore, the analysis of Nilsson-Uggla should be 
extended to the piecewise polytropic balls when we need an accurate interval of the index 
n cor - However, in the following analyses, we simply assume that the physically reasonable 
interval of piecewise polytrope index n cor is the same with the simple polytrope index, 


0 < n cor < 5 . (28b) 

In this paper, we do not aim the thorough analysis of the core-envelope piecewise polytropic 
balls, but we perform a test analysis with the representative model under the extreme-sound- 
speed condition (26) in the parameter region (28). We expect that our model can represent 
some typical behavior of core-envelope piecewise polytropic balls. 

Furthermore, here we introduce one more expectation: Remember that, for the simple 
polytrope (4) as mentioned in section 3, Nilsson-Uggla [13] have revealed the existence of 
divergence of the mass M* and radius 72* at some isolated parameter points in the physically 
interesting parameter region (12). Therefore, for the core-envelope piecewise polytrope (23), 
the same divergence of M* and 72* may occur. However, let us expect that the ratio 3M*/72* 
converges to zero even if M* and 72* of the piecewise polytropic ball diverge, since 3M*/72* 
of simple polytropic ball is expected to converge to zero as shown in figure 1. 


5.2. Problem A: Can a UCOP appear outside the core-envelope piecewise polytropic 
ball? 

The problem which we numerically analyze in this section is whether or not an inequality, 
72* < 3M*, holds for the fluid ball of equation of states (23) under the extreme-sound-speed 
condition (26). Our strategy to calculate the ratio 3M*/72* is as follows: 

Al’: Solve numerically TOV equations (11a) and (lib) for given values of parameters 
(X b ,n CO r), and calculate the ratio, 3M*/72* . 

A2’: Iterate the step Al’ with varying parameters (Xb,n cor ), so as to obtain the ratio, 
3M*/72*, as a function of parameters (Xb,n cor ). 

A3’: Find the maximum value of 3M*/72* as a function of (Xb,n cor )- If the maximum is less 
than unity, we conclude that the inequality, 3M* < 72*, holds for all values of (Xb,n cor ) 
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contours of 0.972 0.974 0.976 



Fig. 6 Contours of 3M*/on Xb-n cor plane for the core-envelope piecewise polytropic 
ball. Numerical calculation is performed in the region, 0.05 < < 1.0 and 0.2 < n < 6. 

Contours of 3M*/R* are plotted at intervals of 0.05 in left panel and at intervals of 0.002 in 
right panel. It is found that the inequality, 3< 1, holds. 

in the parameter region (28), and no UCOP appears outside the core-envelope piecewise 
polytropic ball under the extreme-sound-speed condition (26). 

Our numerical result is shown in figure 6, in which the contours of 3are plotted on 
^tr^cor plane. Because the contour of 3= 1 does not appear, the inequality, 3M* < 
R *, holds for our model of core-envelope piecewise polytropic ball under the extreme-sound- 
speed condition (26) in our interesting parameter region (28). 

Although the numerical analysis shown in figure 6 has been performed under the 
extreme-sound-speed condition (26), our original interest is the case under the subluminal- 
sound-speed condition (25). In order to consider the subluminal case, let us remember the 
discussion for introducing the extreme-sound-speed condition (26). Then, according to our 
previous paper [8], we can expect that no UCOP appears outside the core-envelope piece- 
wise polytropic balls under the subluminal-sound-speed condition (25) in our interesting 
parameter region (28). 

5.3. Problem B: Can a UCOP appear inside the core-envelope piecewise polytropic 
ball? 

The problem which we numerically analyze in this section is whether or not a UCOP can 
exist inside the fluid ball of equation of states (23) under the extreme-sound-speed condi¬ 
tion (26). The calculation and discussion until equation (22) in section 4 is applicable to the 
core-envelope piecewise polytropic ball. Therefore, our task is to check whether or not the 
existence conditions (21) hold. 
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Fig. 7 Plots of curves Ci(E, Eb, ^cor) = 0 and C 2 (X, Xb, ^cor) = 0 in E-Xb plane at poly¬ 
trope index n cor = 3. The region, C 2 < 0, is also denoted. It is recognized that the existence 
conditions of UCOP (21) are not satisfied for this case. 


Note that, for the core-envelope piecewise polytropic ball under the extreme-sound-speed 
condition (26), the quantities C\ and C 2 in equation ( 22 ) can be regarded as functions 
of three arguments, Ci(X, Xb, n cor ) and C 2 (X, Eb, n cor ), since the values of C\ and C 2 are 
calculated by substituting the solution of TOV equations (11a) and (lib). Then, our strategy 
is as follows: 

Bl’: Solve numerically TOV equations (11a) and (lib) for given value of parameters 
(£b,n cor ), and iterate this numerical calculation with varying Xb and fixing n cor at 
a given value. This iteration produces Ci(E, Xb, n cor ) and C^X, Xb, n cor ) as functions 
of (X, Eb) for the given value of n cor . 

B2’: Plot two curves, C\ = (1 and C2 = 0, and identify two regions, C2 > 0 and C2 < 0, in 
X-Xb plane for the given n cor . If the curve C\ = 0 does not intersect with the region 
C*2 > 0, it is concluded that, at the given value of n cor , no UCOP exists inside the 
core-envelope piecewise polytropic ball under the extreme-sound-speed condition (26). 
B3’: Iterate the steps Bl’ and B2’ with varying n cor , and check whether or not the intersection 
of C\ = 0 with C2 > 0 exists at each value of n cor - If the intersection does not appear 
for any value of n cor in our interesting parameter region (28), then we conclude that a 
UCOP can never appear inside the core-envelope piecewise polytropic ball under the 
extreme-sound-speed condition (26). 

By the above strategy, it is found that the quantity Ci(X, Eb, n cor ) remains positive and 
does not become zero {C\ > 0 ) in X-Eb plane for low values of polytrope index, 0 < n cor < 
^cor-low, where n C 0 r _i 0 w ~ 1.5. And, for higher values of polytrope index, n cor -iow < ^cor < 5, 
although the quantity C\ becomes zero at some values of E and Eb in X-Xb plane at every 
value of n cor , the curve C\ = 0 remains in the region C 2 < 0. The example at n cor = 3 is 
shown in figure 7. Hence, we can conclude numerically that no UCOP appears inside the 
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core-envelope piecewise polytropic balls under the extreme-sound-speed condition (26) in 
our interesting parameter region (28). 

Although the above numerical analysis has been performed under the extreme-sound- 
speed condition (26), our original interest is the case under the subluminal-sound-speed 
condition (25). In order to consider the subluminal case, let us refer to the simple polytropic 
balls analysed in section 4, and note that the curve C\ = 0 in £- V c plane for the simple 
polytrope, shown in figure 4, appears in the region of high values of V c . Hence, if the sim¬ 
ilar behavior is expected for the core-envelope piecewise polytropic balls, the curve C i = 0 
in S-Sb plane tends to disappear as the sound speeds Vy, and V c decrease. That is, the 
intersection of Cy = 0 with C *2 < 0 does not appear under not only the extreme-sound-speed 
condition (26) but also the subluminal-sound-speed condition (25). This discussion makes 
us expect that no UCOP appears inside the core-envelope piecewise polytropic balls under 
the subluminal-sound-speed condition (25) in our interesting parameter region (28). 

6. Summary and discussions 

We have investigated whether or not a UCOP can exist in the spacetime of static spherical 
ball of perfect fluid. The equations of states we have considered are the simple polytrope (4) 
and the core-envelope piecewise polytrope (23). By numerical analyses of TOV equations 
and null geodesic equations, our result is as follows: For the simple polytropic balls, we have 
performed the thorough numerical investigation, and concluded that no UCOP can exist 
inside nor outside any simple polytropic ball under the subluminal-sound-speed condition. 
For the core-envelope piecewise polytropic balls, we have numerically investigated the rep¬ 
resentative model under the extreme-sound-speed condition (26), and concluded again that 
no UCOP can exist inside nor outside the core-envelope piecewise polytropic ball under the 
extreme-sound-speed condition. Further, according to our previous paper [8] and section 4 
of this paper, it is expected that the conclusion under the extreme-sound-speed condition is 
also true of the case under the subluminal-sound-speed condition. 

Above conclusions mean that the polytropic balls investigated in this paper cannot be a 
black hole mimicker which possesses UCOP but no black hole horizon. This implies that, if 
the polytrope treated in this paper is a good model of the stellar matter in compact objects, 
the detection of shadow image by optical observation is regarded as a good evidence of black 
hole existence. Note that, to obtain a more definite conclusion for the core-envelope piece- 
wise polytropic ball under the subluminal-sound-speed condition, a more detailed numerical 
research is necessary as discussed in the last two paragraphs of section 5.1. 

Next, let us discuss a by-product of our analysis. In sections 3 and 5.2, the ratio of total 
mass to surface radius of polytropic ball, 3M*/i?*, has been the central issue. As mentioned in 
section 1, the mass-to-radius ratio must be bounded above, 3 M*/R* < 3/2, in order to avoid 
gravitational collapse. Buchdahl [14] decreased the upper bound to, 3M*/i?* < 4/3, with 
assuming non-increasing mass density in outward direction and barotropic equation of states. 
Next, Barraco and Hamity [15] decreased the Buchdahl’s upper bound to, 3M*/i?* < 9/8, by 
adding dominant energy condition to Buchdahl’s assumptions. Furthermore, in our previous 
paper [8], we decreased the Barraco-Hamity’s upper bound to, 3M*/i?* < 1.0909209, by 
replacing the dominant energy condition with the subluminal-sound-speed condition. All 
these upper bounds remained greater than unity, which permits the existence of some black 
hole mimicker. However, as shown in figure 2 of this paper, the upper bound of 3M*/i2* 
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Fig. 8 Sectioned diagram of figure 2 at V c = 1. f(n) describes 3M*/R* of the simple poly¬ 
tropic ball as a function of n at V c = 1. Differential of f(n) is plotted by an approximation, 
d/(n)/dn ~ [f(n + 5) — f(n — S)]/5, where 5 = 10 -3 . 

is decreased to a value lower than unity, by restricting the equation of states to the simple 
polytrope (4) and assuming the subluminal-sound-speed condition (10). Since the upper 
bound is found on the vertical line at V c = 1 in figure 2, a sectioned diagram of figure 2 at 
V c = 1 is useful to read a precise value of the upper bound. It is shown in figure 8, where 

we define f(n) by 3M*/R* as a function of n at V c = 1. From the figure 8, it is concluded 

that the following inequality holds in the physically interesting parameter region (12) for 
the simple polytropic balls under the subluminal-sound-speed condition, 

< 0.844, (29) 

-it* 

where the upper bound is given by parameters V c = 1 and n — 0.78, and the value of upper 
bound is numerically calculated, /(0.78) ~ 0.844. On the other hand, when we restrict the 
equation of states to the core-envelope piecewise polytrope (23) and assuming the extreme- 
sound-speed condition (26), the upper bound of the ratio 3can be roughly read from 
the figure 6, 

—^< 0.976. (30) 

it* 

Finally, let us make a comment on a related topic. A possibility of trapping gravitational 
waves inside stellar objects has been discussed [16, 17]. The potential of gravitational pertur¬ 
bation is analysed in those discussions, while the potential of light propagation such as shown 
in figure 3 is discussed in our analysis. So, when one is interested in a combination of optical 
observation and gravitational wave observation in the search of black holes and gravitational 
waves, the existence/non-existence conditions of UCOP and gravitational-wave-trapping 
may become an interesting issue. 
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A. On numerical treatment of TOV equations 

Right-hand sides of TOV equations (11a) and (lib) are indeterminate form at center 
because of the conditions, M —>• 0 and R —>• 0 as £ —>• 1. Therefore, in solving TOV equations 
numerically, we have made use of perturbative solutions near the center. 

In order to consider a perturbation near center, we regard the radius R as an independent 
variable, and the mass density as a function of radius, £(P). TOV equations (11a) and (lib) 
are rearranged to 


d M(R) 
d R 

d P(R) 
d R 


4ttP 2 £(P) 

[£(P) + P(P)] [ M(R) + 4itR 3 P(R) ] 
R [R - 2M(R) ] 


(Al) 


For a sufficiently small radius R <C 1, we introduce perturbations, 


M(R) = M(\)R + M( 2 )R 2 + M( 3 )R 3 H- 

P(R) = P C + P(i)R + P{2)R 2 + P(3)R 3 + • • • (A2) 

£(P) = 1 + £(i)P + £(2 )R 2 + S( 3 )f ? 3 + • • ■ , 


where conditions M(R = 0) = 0, £(P = 0) = 1 and P(R = 0) = P c are included. Substitut¬ 
ing (A2) into (Al), we obtain Mnj = 0, M 2 = 0, Po) = 0 and remaining parts, 

M(R) = -nR 3 + vrS (1) P 4 + • • • 

3 

P(R) = P c - |tt( 1 + 3P c ) (1 + P C )P 2 - ^(7 + 15P C )S (1) P 3 + • • • ( A3 ) 

o y 

S(P) = 1 + + S(2)P^ + P‘(3)R 3 + • ■ • 1 


where the central pressure P c and coefficients (n = 1,2, 3, • • •) are determined by 

concrete form of equation of states. 

Substitute these perturbative expansions into the equation of states (4), we obtain 


^ n ^ n 27t( 1 + P c ) (1 + 3P C ) 

£ (1) = 0 , E(2) = -rp- 


(A4) 


Hence, denoting a small radius by < 1, the mass density S 5 and mass Ms at R = Rs 
are approximately given by £5 = 1 + S( 2 )P 2 and Ms = (47r/3)P 3 . If the mass density near 
center S5 are given, then the others are determined by 


Rs 


1 — V,5 

|S(2)| 


M S 


4/ 1 — Sa 3/2 
3 V |S( 2 )| 7 


(A5) 


In numerical calculation, we have solved TOV equations (11a) and (lib) for interval 0 < S < 
£5 with initial condition (A5). Also, we have checked the convergence of numerical solutions 
with varying £^. All results in this paper are obtained using £5 = 1 — 10 -4 . 
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